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Abstract 

Quantum Chromodynamics on a lattice with Wilson fermions and a chirally twisted 



Qh! mass term is considered in the framework of chiral perturbation theory. For two and 

three numbers of quark flavours, respectively, with non-degenerate quark masses the pseu- 
doscalar meson masses and decay constants are calculated in next-to-leading order includ- 

• i— i . 

ing lattice effects quadratic in the lattice spacing a. 

a' 

When quantum field theories are simulated numerically, the introduction of a finite space-time 
lattice brings with it two standard problems. In order to obtain results relevant to the physics 
in the infinite continuum the effects of the finite lattice spacing a and of the finite volume of the 
system have to be controlled. The resulting demands on the lattice are limited by the available 
computer resources. Simulations of Quantum Chromodynamics are facing a third difficulty. 
The common algorithms slow down severely when the masses of the light quarks are decreased 
towards their physical values. 

Recent investigations [H El E| suggest that lattice QCD with Wilson fermions and a chirally 
twisted quark mass term jH [H] offers an attractive framework for numerical simulations with 
small up- and down-quark masses. In particular, for the case of full twist lattice artifacts of 
order a are removed automatically For a review of twisted mass lattice QCD see [2j. 
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An important tool for QCD is chiral perturbation theory [HI 03 EDI- For sufficiently small 
quark masses it describes the dependence of mesonic quantities like masses and decay constants 
on the quark masses in terms of expansions in powers of quark masses (modified by logarithms). 
These expansions contain low energy parameters, the Gasser-Leutwyler coefficients, which have 
to be determined by methods outside chiral perturbation theory. 

Chiral perturbation theory is particularly suited for lattice QCD. In its range of applica- 
bility, it can be used to extrapolate the results of Monte Carlo simulations into the region of 
small quark masses. On the other hand, numerical simulations of lattice QCD can provide 
information about the Gasser-Leutwyler coefficients, see [TT] for a review. 

The adaption of chiral perturbation theory to lattice QCD |T21 H31 UU EEH1 dEj brings the 
lattice spacing a in as a second expansion parameter. In the effective chiral Lagrangian lat- 
tice artifacts are represented by additional terms proportional to powers of a [TJ]. Physical 
quantities appear in double expansions in quark masses and in a. 

Chiral perturbation theory has been extended to the case of twisted mass lattice QCD in 
[I3IIE1- The resulting expansions contain the twist angle as an additional parameter. Results 
for various quantities have been calculated in next-to-leading order of chiral perturbation theory 
including lattice effects of order a 2 [13 [HI CHI EDI EH E21 EHl EU . 

In most of the work the case of degenerate quark masses is considered. Current numerical 
work on twisted mass lattice QCD with non-degenerate quark masses [2E| motivates to imple- 
ment this situation in chiral perturbation theory. This is the purpose of the work presented 
here. 

We consider lattice QCD with Nf = 2 or 3 quark flavours and a twisted mass term. In 
the two-flavour case we consider split quark masses (m u ^ m^). In the case of three flavours 
we restrict ourselves in this letter to equal up- and down-quark masses and a heavier strange 
quark mass (m u = m d ^ m s ). The masses and decay constants of the pseudoscalar mesons are 
calculated in next-to-leading (one-loop) order, including lattice terms of order a 2 . 

Chiral perturbation theory for pseudoscalar meson masses in untwisted lattice QCD with 
Nf = 3 flavours of Wilson quarks has been considered in [26J. These authors use a power 
counting scheme different from the one of [TJ] [TH H3 [THJ and a different parameterization 
of the action, which hampers a comparison of our calculations. The results which could be 
compared in the limit of vanishing twist angle do agree. The case of twisted mass lattice QCD 
with 2+2 flavours has been studied recently by means of chiral perturbation theory in [22j. The 
implementation of mass non-degenerate quarks in twisted mass chiral perturbation theory for 
two flavours has been considered by Walker-Loud and Wu [2H]. They extend the formalism to 
the baryon sector and calculate nucleon and delta masses. 

In lattice QCD a chiral twist of the quark mass term is generally introduced by transforming 
the quark mass matrix M into 

M{uj) = e iuJbTb ^ Me iUiT ' 75 , (1) 



2 



where Tf, are the generators of SU(iVy) and uj denotes the set of group parameters Ub- In 
contrast to the situation in the continuum, on the lattice the chiral transformation of the mass 
matrix has a measurable effect on observables since in lattice QCD with Wilson fermions chiral 
symmetry is broken by lattice artifacts. 

In chiral perturbation theory the low energy regime of QCD is described by an effective 
Lagrangian for the pseudo-Goldstone bosons of spontaneously broken chiral symmetry. For 
Nf = 3 they form the pseudoscalar meson octet including the pions, kaons and r/. 

The transformation properties of the Goldstone bosons in QCD lead to collect the dynamical 
variables in the SU (Nf) matrix valued field 

U(x) =exp (^2yT b M^ , (2) 

where are the components of the meson fields and F is a low energy parameter. Treating 
quark masses and lattice effects perturbatively in a spurion analysis with respect to the required 
symmetry properties, the leading order effective chiral Lagrangian for lattice QCD is given 

by nam 

£ LO = £2 = Tr {dJJ d^) - S Tr ( x rf + U X ^) -^Tr(ptf + Up^) , (3) 

consisting of the chiral symmetric kinetic term and two chiral symmetry breaking terms. The 
first of them contains the quark mass term, given by 

X = 2B M, (4) 

and in the second one lattice artifacts are represented by 

p = 2W at, (5) 

where B and Wq are additional parameters. 

Note that we use the power-counting scheme of [HUE], accounting for terms of order {m, a} 
in leading order and {m 2 ,ma,a 2 } in next-to-leading order. In order to define an applicable 
scheme we always assume the lattice effects to be smaller: W a <C B m q . 

Introducing a chiral twist, the mass matrix in the effective Lagrangian becomes 

XH = 2B e - iuj » n Me- [uJbTb , (6) 

in the so-called twisted basis. The (^-dependence of the mass term can be removed by changing 
to the physical basis, i.e. by a chiral transformation of the configuration matrix 

U — > U' = e iuJbTb Ue iuJbTb (7) 

and a corresponding chiral twist of the lattice term according to 

p — > p{oo) =e iu;bTb pe iu,bTb . (8) 
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Omitting the prime on the meson field and suppressing the explicit ^-dependence of the lattice 
term, the effective Lagrangian looks like the one of untwisted lattice QCD. However, the intro- 
duction of a chirally twisted mass term leads to a shift of the minimum 7^ of the effective 
Lagrangian of order a. Chiral perturbation theory, which is based on an expansion around 
this minimum, has to be adapted accordingly. We treat the nontrivial minimum perturbatively 
according to the applied scheme. Denoting the minimum 

U = exp\2yT b &j , (9) 

a substitution 

u — ► u 1/2 u u 1/2 (10) 

is performed before expanding the Lagrangian. The tree-level and one-loop contributions of 
the resulting C 2 have the same form as those of continuum chiral perturbation theory, but with 
different leading order meson masses which include lattice artifacts, see below. 

The next-to-leading order part of the effective Lagrangian in the present power counting 
scheme is given by [HH Ej 

£4 = - U [Tr (d^Ud^)] 2 - L 2 [Tr {d^Ud^)] 2 - L 3 Tr {[d^Ud^f) 

+ L 4 Tr (d^Ud^) Tr ( X ] U + U ] x) + W 4 Tr (d^Ud^) Tr (p f [7 + Wp) 
+ L 5 Tr (d,Ud,U^U X ^ + XU ] )) + W 5 Tr (d,Ud,U^Up^ + p W)) 

- L 6 [Tr ( X ] U + U ] x)] 2 - W 6 Tr ( X ] U + U ] X ) Tr (p'U + W p) - W' & [Tr (p*U + U ] p)] 2 

- L 7 [Tr ( X ] U - U ] x)] 2 - W 7 Tr ( X ] U - W X ) Tr (p ] U - Wp) - W 7 [Tr (p*U - Wp)] 2 

- L 8 Tr {x ] U X ] U + rf X U*x) ~ W 8 Tr (p^U^U + U^p rf X ) 

-W^Ti{p ] Up ] U + U ] pU^p) . (11) 

Next-to-leading order terms are given by tree-level contributions from £4 and one-loop 
contributions from C 2 . The latter produce divergencies, which are treated by renormalizing the 
low energy coefficients L i? Wi and W[. 

First we consider Nf = 2 flavours with non-degenerate quark masses. Following (29|, to 
ensure positivity of the fermion determinant we implement the quark mass splitting and the 
chiral twist with different generators T&. Because T 3 = r 3 /2 is usually chosen for the twist, we 
use Ti = Ti/2 for the mass difference. Defining 

Xo = B (m d + m u ) and Xi = B (m d - m u ), (12) 

we have 

x = Xoi + xin (13) 

and 

p = p t + \p3T3 , (14) 
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where po = 2Wqcl cos uo and p% = 2Woasinu;. Due to the properties of the Pauli matrices only 
the L 7 and L 8 -terms lead to contributions involving the mass splitting \i- 
For the minimum of the effective Lagrangian we obtain 







(15) 



and 



Xo 



1 - t^(2L 86 - W 8e ) X o 



+ 0(a 2 ) , 



where terms of order a 2 can be neglected at this order of the calculation, and we define 

Lij = Li + NfLj . 



(16) 



(17) 



From a one-loop calculation of the renormalized propagator we obtain the pion masses in 
next-to-leading order. Owing to the unequal quark masses and the chiral twist the three masses 
are generally different. For the component <p2 = ^2 we get 



M 2 2 =ml + 



+ 



32tt 2 F 2 V A 2 



F 2 



+ 2U m {xl ~ Pi 



where m n denotes the leading order pion mass, 



L r 5A Xo{xo + Po) ~ W&(XoPo + Po + p\ 
2^ 8 r 6 (xoPo + Pi) + 2W* p 2 



2 1 P% 

r m 7T = Xo + Po + 



2Xo 



(18) 



(19) 



A is the renormalization scale, usually taken to be 47rF , and the superscript r indicates the 
renormalized L, W and W coefficients. 

The masses of the other two pions can be read of from the mass splittings 
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M 2 — M 2 





16 

^2 



(20) 
(21) 



Setting Xi = 0, this is consistent with the results for degenerate quark masses [13 [201 - Moreover 
we observe O (a) -improvement at maximum twist uj = ir/2, where po = 0. 

Note that with our choice of SU(2) generators, ~K\ is the neutral pion, and the charged pions 
are linear combinations of n 2 and 7r 3 . Due to the twist, however, the mass eigenstates are not 
the usual charged pions but the components ir 2 and tt^. 

The pion decay constant F n is given by the matrix element 



(0\A^ a \4> b (p)) = iFjPpSab, 



(22) 
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where A^ a is the axial current. In next-to-leading order one can determine the decay constant 
including order a terms. We get 

,2 



Fn 1 

F 16tt 2 F q 2 



m l 111 ( ) + ^( L 54X0 + Po) 



(23) 



without any contributions of the quark mass difference Xi an d coinciding with the result of 

Now we turn to the consideration of Nf = 3 flavours with degenerate up- and down-quarks 
and a heavier strange-quark. The generators of SU(3) are given by the Gell-Mann matrices 
Tb = Aft/2. Defining Xs = 2B m Sl the quark mass matrix is given by 



X = \( 2 Xo + Xs)t ~ ^(^* _ *o) A 8 • 



(24) 



For Nf = 3 it is generally possible to introduce a twist with two twist angles involving all 
flavours. In view of the fact that the motivation for twist is to eliminate infrared problems 
stemming from the light quark masses, and in order to ensure a positive quark determinant 
[30J we implement a chiral twist by a single angle u in the isospin sector of the lightest quarks, 
leading to 



/po + ips \ 
po - i P3 



V 











1 



(2po+|p|)l 



(|p| - p )A 8 + 1P3A3 



(25) 



with \p\ = 2Wo a. 

The minimum of the potential is shifted again in T3-direction to 



7T 3 



£q — 

Xo 



1 - t^{(2L 8 - W 8 )xo + (2L 6 - W 6 )(2 Xo + Xs)} 



At leading order the pseudoscalar meson masses are 



2 P% 
m; = Xo + Po + 



2Xo ' 



(Xs + Xo) + (\p\ + Po) + 



A 

2xo 



(2x, + Xo) + (2|p|+p ) + 



2xo 



0{a 2 



(26) 



(27) 
(28) 
(29) 



Including the tree-level and one-loop terms at next-to-leading order we obtain the pion 
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masses 



M, 



Tl 2 



2 m TT 



3m^ In 



A 2 



- m v In 



to; 

7 

A 2 



2\ n 



+ 



(4LJ + 8L£ - 2L£ - ALDxl + (4L£ - 2L r 4 ) XoXs 



+ (-2LI - AL\ + mi + mi - 2WI - AWDxm + (2WI - 2Wl)xM 
+ (-2LI + 2Wl) Xs p 

+ (-2W 5 r - 4W^ 4 r + AW£ + 8^ r )p' + (-2WI + 4W^ r )p |p| 

+ (-4LJ - 8LI + AWl + 8WI - 2W r h - AWl)p\ + (-2LJ - + 2WJ)^g] . (30) 

This expression shows that the presence of the strange quark, whose mass term is not subject 
to any chiral twist here, foils the C(a)-improvement at maximal twist, as expected. 
The mass splitting, given by 



M; 



7T1 2 



Ml 



^0 



is formally equivalent to that obtained for Nf 
coefficients into account. 
For the kaons we get 



Wl 6 + W£)pl , (31) 
2, taking the different Gasser-Leutwyler 



M 2 K = m\ + 



m 



K 



A 2 



+ 



+ 



[2L r 8 + 8L r 6 -LI- ALDxl + (AL r 8 + 16L£ - 2L\ - QL\) XoXs 



+ 
+ 



2L; + 4Ll-Ll-2Ll) X 



+ 


1 TT 

y L h 


AL\ + 2WI + 8WI 


-wi 


- AWl) Xo p 


+ 


y L 5 


AL\ + 2WI + 6WI 


-wi 


-2Wl) X ,\p\ 


+ 


y L <o 


2L\ + 2W'l + QWl 


-wi 


- 4Wl) XsPo 


+ 


-Ll- 


2L\ + 2WI + AWl 


-wi 


- 2Wl) Xs \p\ 


+ 


[-2WI 


- QWl + AW' s r + 120o(po + |p|) 



— 4Lg - 8L r 6 + AWl + m l - 2W 5 - m l + aw q)pI 
-2U 8 - 6V 6 + LI + 2WI + 6WI - WI - 4W 4 > 2 .g 



(32) 
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and for the 77-mass 



Ml 



2, m v 

+ 48tt 2 F 2 



3m|- In 



m 



A 



A 2 



- 2m z v In I 



/ 9 

/ mi 



V A 2 



+ 

+ 

+ 
+ 
+ 
+ 
+ 
+ 
+ 



mi 



96tt 2 F 2 



A 2 



A' 



A 2 



9F 2 



[12U 8 + 24L? + 24L r 6 - 2L\ - \2L^xl 



-A8L r 7 + 60L r 6 - 8L r 5 - 30Ll) X oXs + (2AL r 8 + 2AL r 7 + 2AU 6 - 8L; - 12L\) X 

-2L; - \2L\ + Y2WI + 24Wf + 24W 6 r - 2W r h - 12W 4 r )xoPo 

-4L£ - 24L^ - 24W 7 r + - 4W 5 r - 6W 4 r )x |p| 

-4L T 5 - 6^ - 24W 7 r + 30W£ - 4W 5 r - 24W 4 r )x s p 

-8L£ - 12L£ + 24W 8 r + 24W 7 r + 24W 6 r - 8W; - \2Wl) X s\p\ 

-10W 5 r - 24WI + 36W£" + 48Wf + 48^ r )Po 

-8W; - 30WX - 48Wf + 60^6 r )Po|p| + (-24L^ - 2AL r 7 - 24L£ 



+ 24W 8 r + 2AW r 7 + 24W 6 r - 10W 5 r - 24W£ + 12W^ r + + 24W 6 ' r )p§ 

+ (24L; - 30L£ + - 24W 7 r + 30W 6 r - 4W 5 r - 24W 4 r )p 2 ^ . 



(33) 



In addition to the masses we determine the decay constants for the pseudoscalar meson 
octet. The results are for the pions 



for the kaons 



F 



1 - 



32tt 2 F 2 



27n In 



7774 
A 2 



+ 777 K In 



77). 



A 



A 2 



+ --2 [Llxo + LKxs + 2 X o) + W^po + W£(\p\ + 2p )] , 
^0 



(34) 



Fk 
F 



1 - 



128tt 2 F 2 



771 In 



rn 



+ -^[Ll(xo + Xs) + 2n 4 (xs + 2 Xo ) 
^0 



and for the 77 meson 
3 

= 1 

^0 



^ = 1- 
F n 32tt 2 F 2 



777 ^ In 



rn 



K 



A 2 



/// 



A 



A 2 



+ m rj In 



777^ 
A 2 



2\ 1 



Wa\p\+p ) + 2Wl(\p\+2p )}, 



(35) 



+ — T2 \L\{2 Xs + Xo ) + SLUxs + 2 X o) + W£(2\p\ + p ) + 3W^\p\ + 2p )] . (36) 



In the limit of equal quark masses for all flavours, x 
expressions coincide to 

,2 



Xo, and vanishing twist these 



Y = 1 

^0 



N 



f 



32tt 2 F 2 
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rn; 



4 

+ ~B2 t L 54X0 + W54A)] > 
^0 



(37) 



recovering the result of [H]. 

As a check, in the continuum limit, p = 0, our formulae for Nf = 3 go over into those of 

M 

The expressions given here can be used in the analysis of numerical results from unquenched 
simulations of twisted mass lattice QCD with two or three flavours. Fitting the data for masses 
and decay constants with the chiral perturbation formulae will allow to extract numerical values 
for the low energy coefficients |3| 125]. 

We thank I. Montvay, F. Farchioni and P. Hofmann for discussions. 
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